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CLT Abstract 
<D ' 

We calculate the S-wave charmonium spectroscopy using the Hamiltonian with the non- 
CN ■ relativistic QCD (NRQCD) potential. The logarithmic factor In fir, appearing in the next- 

to- leading order QCD loop corrections to the potential, is expanded about r = 1/fi, where 
fi corresponds to the typical charmonium scale. The resulting potential characterized by the 
Coulombic and linear components is consistent with the form of the Cornell potential. We ob- 
tain x 2 fitting results for the masses of the S'-wave charmonium states, r; c (l 1 1 S , o), J/?p(l 3 Si), 
rj c (2 1 So), and ip(2 3 Si) in remarkable accordance with data. Our results successfully account for 
the hyperfme splitting for the IS state as well as for the 2S state. We further use the three best 
fit parameters: the charm quark mass m c , coupling constant a s and the corresponding scale \i 
^ . to predict the 5-wave mass spectrum with n < 6. The hints for results are discussed. 
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I. INTRODUCTION 



The heavy quarkonium (like charmonium, bottomonium, etc.) is a system where we can 
study the low-energy QCD in a systematic way. The heavy quarkonium satisfies the following 
hierarchy scales 

M»p~l/r~MO>£~ Mv 2 , 

where M is the heavy-quark mass, p the momentum transfer, E the binding energy, r the 
typical distance between the quarks, and v the typical heavy quark velocity. It can be estimated 
that v/c ~ 0.6 for charmonium and 1/3 for bottomonium Q], Taking into account the 
above properties, one can introduce the nonrelativistic effective field theory which realizes a 
factorization at the Lagrangian between the high energy effects and low energy contributions. 
By integrating out the hard parts we can obtain the non-relativistic QCD (NRQCD) which is 
expanded in p/M and E/M [l|, 0, Q|. Taking into account the fact that for the charmonia and 
bottomonia the typical scale \i ~ p associated to the inverse size of the system 1/r is satisfied by 
the relation 1/r) ^> E ~ Aqcd, the NRQCD can be further expanded in Er and leads to 
an effective field theory which is the so-called potential non-relativistic QCD (pNRQCD) 0-0]. 

Since the QCD has been widely accepted as a fundamental theory for the strong interactions, 
we have no doubt about that it should describe the spectroscopy of the heavy quarkonium. Nev- 
ertheless, in practice, it is still not so successful. On the other hand, the phenomenological quark 
models which mimic the QCD features seem to offer better results for the heavy quarkonium 



systems. The form of Cornell potential used in the phenomenological quark model IsKioi], where 
the potential is made of the Coulombic and linear parts, has been confirmed to be valid by the 
lattice calculation [ll|) 0] • 

The cc charmonium states are usually denoted by the symbol n 2s+l Lj with n and s being 
the principal and total spin quantum numbers, respectively. The r/ c (2 1 S'o) was first measured in 
B decays by Belle in 2002. Most potential model calculations predicted a much low value for 
its mass compared with the data. On the other hand, so far the hyperfine splitting for the 2S 
state (?n^(2 3 Si) ~~ m vc(2 1 S ) = 49 ± 4 MeV) as well as for the IS* state (mj/^nag^ ~ m "n c {i 1 So) = 
116.6 ± 1.2 MeV) cannot be simultaneously calculated well compared with the data. Although, 



the phenomenologically potential models [9j, |l0j, Il3l4l6| can offer quite intuitive picture for cc 



charmonium spectroscopy, the deviation between the theoretical calculations and data is quite 
large, so that we cannot make further predictions accordingly. Even the old ■0(4415), which 
was conventionally assigned as the V(4 3 5i) quantum state, has been argued that it should be 
V>(5 3 Si) 01! 

Recently, a number of interesting charmonium-like states, that are above the DD open- 
charm mass thresholds and named collectively as "X, Y,Z" mesons, have been found (see the 
discussions, e.g., in Refs. 0, Q]). So far, most of them are at odds with expectations of the cc 
states which had been predicted by conventional charmonium models. One may then suggest 
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that some of the "X, Y,Z" particles are exotic. Existence of exotic states such as glueballs, 
hybrid mesons (the bound states of ccg), molecules, and tetraquark mesons (the bound states of 
cccc), which go beyond the description of the naive quark model, can offer the direct evidence 
concerning the confinement property of QCD. Nevertheless, so far, none of the exotic states can 
be well established. On the other hand, it is still difficult to assign any new observable in the 
conventional cc charmonium mass spectrum with the definite quantum state; many suggestions 
can be found in the literature. For instance, it was suggested that the X(3940) may be the singlet 
state r] c (3 So). However the corresponding triplet state 3 3 Si is ^(4040) with mass 4039±1 MeV, 
so that the assignment for V(3940) implies a larger singlet-triplet mass splitting (~ 100 MeV) 
for radial number n = 3 than that (~ 50 MeV) for n = 2, which is one of the problems. 

Motivated by the above reasons, in the present study we will try to obtain an evaluation for 
the S-wave charmonium spectroscopy starting with the pNRQCD Hamiltonian, instead of the 
phenomenologically potential model. We hope to clarify some ambiguities between observables 
and theoretical calculations. We use the QCD potential Vqcd(O) which can be obtained from 
matching NRQCD to pNRQCD 0]. To solve the mass spectrum, we expand In fir, resulting from 
the QCD loop corrections to the potential, about r = 1/^, where [i corresponds to the typical 
charmonium scale of order m c v. The benefit of the expansion is that our resulting potential 
exhibits the form of the Cornell potential which was confirmed by the lattice calculation [l^ . 
In our study, we have three parameters, the charm quark mass, a s and /i, which are related 
to the determination of charmonium masses. Using the modified Pade approximation [ig| . I20I ] 
which is a numerical technique, we perform the best x 2 -fit for the current data of masses of the 
well-established 5-wave charmonium states, r/ c (l 1 5o), J/ip(l 3 Si), ?7 c (2 1 5o), and ip(2 3 Si), and 
then use the fitted parameters, m c , a s , and \i to further predict the full S- wave mass spectrum. 
^(2 3 Si) is usually denoted as the state ^(3683). As the fact that the total angular momentum 
J is a conservation quantum number, the orbital angular momentum L is actually not a good 
quantum number for the cc charmonium states. Therefore t/>(3683) and ^(3777) could be the 
mixtures of 2 3 Si and l 3 Di states 21[. Our result shows that the minimum % 2 is consistent with 
zero, which hints that the S — D mixing may be negligible. 

The Pade approximation is to approximate a function fix), which is expanded in a Taylor 
series up to order k, by the ratio of two polynomials, one of order M in the numerator, and 



another of order N in the denominator, with M + N = k [22J, |23|]. This ratio is called the Pade 
approximant of f(x). The technique of the Pade approximation has the following advantages. 
First, it can accelerate the convergence of the usual Taylor expansion for a given function. 
Second, even for x going beyond the radius of convergence of the Taylor's series of a given 
function f(x), its Pade approximant could well approximate the original function, i.e., physically 
it can be applied to the non-perturbative region. This method thus has been exploited in 



statistical physics, hadron phenomenology, quantum field theory 244291]. and recently in finding 
the solutions of general relativity [30(]. 
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The Pade interpolation method, which is called the modified Pade approximation here, was 



first proposed in Refs. 19|, |20[] to explore physics in the non-perturbative region. In this modified 
approach, a single Pade approximant is obtained by interpolating the weak and strong behav- 
iors. We adopt this approach to study the charmonium spectroscopy. The QCD Hamiltonian is 
redefined as H(j3) = Hq + (3Hl, where He involves the Coulomb-like potential and Hl contains 
the linear potential. We introduce the parameter b to separate the kinetic energy term into two 
parts and then lump into He and Hl separately. As /3 = 1, we have H(l) = H, the physical 
Hamiltonian. We consider two limits, /3 > 1 and /3 <C 1, to perform the perturbation calcula- 
tion. After obtaining the results in the two limits, we can then get the physical eigenenergies 
corresponding to physical Hamiltonian H using the Pade interpolation. (See Sec. IIIIBI for the 
details.) In performing the fit, we also put the constraint on b, so that the numerical error in 
the approach due to the choice of b is small enough (< 2%). In general, when the radial number 
n < 6, the error is less than 7% for 0.1 < b < 0.6. The detailed discussion for numerical errors 
will be presented in Sec IIV1 

The remaining of this paper is organized as follows. Together with an example, we will give 
a brief introduction to the methods of the conventional and modified Pade approximations in 
Sec. HU We formulate the modified Pade approximant for charmonium masses in Sec. IIIH In 
Sec. HVl the prediction for the 5-wave charmonium mass spectrum, together with the best fit 
parameters, m c , a s , and fi, are given by minimizing x 2 fit- Sec. |V]is our summary. 

II. THE PADE APPROXIMATION 

A. The conventional Pade approximation 

The Pade approximant f[M/N](x) of degree (M,N), developed by H. Pade, is an approxi- 
mation of a given function f(x) as a ratio of two power series: 

f[M/N](x) = ^4 (2.1) 
Qn(x) 



PO + PlX + P2X + P3X -\ h PMX 



(2.2) 



1 + qix + q2X 2 + q3X 3 H h pnx n 

where Pm(x) and Qn(x) are polynomials of degrees M and N, respectively. Assume that f{x) 
is analytic around x = and has the Taylor expansion (or called the Maclaurin expansion) form: 

oo 

f(x) = Y,a t x l . (2.3) 

i=0 

Setting /( n )(0) = f^[M/N}(0) with n = 0, 1, . . . , M + N, one has 

M \ / N \ 

E ~ = °(^ +M+1 )' (2-4) 
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which can lead to ./V + M + 1 linear equations: 

ao ~ Po = 
qia + a\ - p\ = 



and 



QM a N—M + Qm-io-n-m+i H 1- oat - pat = (2.5) 

QMO-N-M+l + <lM-ia>N-M+2 + • • • + O'lOAT + OftT+l = 
QM&N—M+2 + QM—l a N—M+3 + ' ' ' + <Zl a iV+l + = 



QMOTV + ^Af-lfflAf+l + • • • + giOjV+M-l + O-N+M = (2.6) 

From the above independent equations, the N + M + 1 coefficients, pi and qi, can thus be 
determined. 

For a given analytic function, its Pade approximant of degree (M, N) often gives much better 
approximation than truncating its Taylor series of degree M + N, and, moreover, the former 
may still work when the latter does not converge. Physically, this implies that not only the 
perturbative results can be further improved, but also it becomes possible to obtain a good 
estimate for nonperturbative phenomenologies. 



B. The modified Pade approximation 

In a practical calculation, we may not know well the full Taylor expansion of a given physical 
quantity at the specific point, e.g., x = 0, but just have its series up to a typical order. Following 
the idea by Leung and Murakowski [lsj, the Pade approximant of the function can be further 
improved if we know the truncated Taylor series of this function at the other analytic point. 
Here we would like to define the modified Pade approximant for a given function as follows. For 
an analytic function /(a?) in the considered range of variable x, if we know its truncated Taylor 
series of degrees r and s respectively at x = and x = xq / 0, 



/Taylor (x) = ^ + 0(x r+1 ), 

i=0 

s 

/Taylor^) = ]T b t { X - XqY + 0( X S+1 ), (2.7) 
i=0 

in analogy to the relation given in Eq. (|2.4|) . we can obtain r + s + 2(= M + N + 1) independent 
equations to determine the modified Pade approximant f( r ' s ^[M/N](x): 

J [M/lV m 1 + qiX + q2X 2 +q3X 3 + ... +pNX N » V-X) 
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which may provide an accurate estimation for the original function in the entire range between 
the two expanding points. Here we take the function f(x) = ln(x + 1) as an example to 
illustrate the points. Expanding about the origin, which is equivalent to modeling the physically 
perturbation, the Taylor series of this function reads 



/ Tay io I .(x) = 5:(-ir 1 -, (2.£ 

1=1 

which converges only for — 1 < x < 1. We can thus get the conventional Pade approximants, 

/[2/2](x) 



T 2 

7C +X + l 



fwm = h^ x ' +x - (2.io) 

4 + 1 



On the other hand, we perform the Taylor expansion for the function at a large value of x, e.g. 
x = 6, which is equivalent to the case of modeling the extremely strong coupling, reads 

/ Taylor (x) = log(7) + ^ - l(x - 6) 2 + ^ - ^ + 0((x- 6f) . (2.11) 

From the Taylor series results given in Eqs. (12. 9h and (12. lip , we obtain the modified Pade 
approximants: 

3(-6+71og(7)) + X 

(l32-71o g (7)- ^^mi),3 1 / 12 71 og(7) N , 

I /I l ^ _0/l_L7 1r.rr('7l _ 0/1 _1_7 1 ^nV7"l ) X T X 



f (2,l)ro/H^\ _ 3024 T ^ I -24+71og(7) -24+71og(7) 

J W^JW - (-12-71og(7))x . 1 

4(-24+71og(7)) + 1 

In Fig. [Hwe plot the exact function f(x), f[2/2](x), f[3/l](x), and / (2,1) [3/1] (a?) , together with 
the result of truncated /Taylor (x) expanding at x = up to 0(x 5 ). More detailed numerical 
results are listed in Table HI The Taylor polynomial (corresponding to the perturbation) is 
valid only for — 1 < x < 1. The results of the Pade approximations can offer reliable estimates 
extending to s > 1, corresponding to the non-perturbative region, and the accuracy can be 
enhanced if one increases the degree(s), M and/or N, of the Pade approximant. The modified 
Pade approximant /( 2,1 )[3/l](:c), which interpolates the results of the two different expanding 
points, differs from the exact value by no more than 1% error between x = and 6. Even for 
/( 1 ' 1 )[2/l](x), the error is still less than 3%. Extending to /^ 2 ' 2 ^[3/2](x), the error becomes less 
than 0.1%. 
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1. 2. 3. 4. 5. 6. 7. 

X 

FIG. 1: The graphs of the exact f(x) = ln(l + x) (solid curve), its Pade approximants, f[2/2](x) 
(long-dashed curve) and /[3/l](x) (dot-dashed curve), modified Pade approximant /( 2,1 )[3/l](x) 
(dotted curve), and Taylor's polynomial of degree 5 (short-dashed curve). 

TABLE I: Comparison of the exact solution f(x) = ln(l + x) and its Pade approximants as 
well as the Taylor polynomial, expanding at x = and up to 0(x 5 ), where the modified Pade 
approximants corresponding to the Taylor expansions both at x = and 6. 



X 


f(x) (Exact) 


/[2/2](x) 




/W)[2/l](x) 


/(^ [3/1] (X) 


/Taylor to 0(x 5 ) 


0.5 


0.4055 


0.4054 


0.4053 


0.4146 


0.4040 


0.4073 


1 


0.6931 


0.6923 


0.6905 


0.7116 


0.6878 


0.7833 


2 


1.0986 


1.0909 


1.0667 


1.1209 


1.0872 


5.07 


3 


1.3863 


1.3636 


1.2692 


1.4021 


1.3751 


35.85 


4 


1.6094 


1.5652 


1.3333 


1.6172 


1.6026 


158.13 


5 


1.7918 


1.7213 


1.2719 


1.7938 


1.7896 


502.92 


6 


1.9459 


1.8462 


1.0909 


1.9459 


1.9459 


1291.2 



III. FORMULATIONS OF HEAVY QUARKONIUM MASSES 
A. The Hamiltonian for cc bound states 

The Hamiltonian for the cc system expanding both in a s and in l/m c , determined from the 
perturbative QCD, is 

H = 2m c + # (0) + V ret (3.1) 
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where m c is the mass of the charm quark at the scale \i. is 

# (0) = — + V 8 , 
including the kinetic energy and static potential up to 0{a 2 s ) 

V - Cpdls _ Cf/3qQ^ In rjx 
r 2tt r 

where P is the momentum of the charm quark, a s the strong coupling constant, and 



(3.2) 
(3.3) 



C F 



, , ai+7 E /3 /2 
I H a. 



7T 



4 „ 2n f 



with 



32 



(3.4) 



(3.5) 



ai = 1.75, rij = 3 for the charmonium 

and "fE — 0.577216 being the Euler constant. We consider the spin-dependently and spin- 
independently relativistic corrections V re i up to 0(a s /m 2 ) and 0(l/m^), respectively, 

Vrel = V LS + VT + V hf + V rel , K . (3.6) 

Here Vls, Vt and Vhf are the spin-orbit, tensor and hyperfine (i.e. spin-spin) interactions, 
respectively, and V re i k is the relativistically kinetic correction, which are given by 

3C FCt s 



V LS 
V t 
Vhf 



V 



reLK 



2 m 2 r 3 

Am 2 r 3 
8itCfO! s 
3m 2 
1 P 4 
4 mf. 



L-S. 



S\2 

Si • s 2 <y 3 (f 



vi, c o(^- Si)(S 2 • f) ^ 
with S12 = 3 ^ Si • 62 , 



(3.7) 

(3.8) 
(3.9) 
(3.10) 



To solve the charmonium spectroscopy, we approximate the Taylor expansion of In r at r 



In fj,r 



1 ( 1 
— r 

2 V H 



n 2 + o 



i 1 



3 loo 
-- + 2r/i- -r /i 



(3.11) 



(3.12) 



with truncated series of degree 2. For the charmonium, // is the typical charmonium scale of 
order 1/r ~ m c v, where v is the velocity of the charm quark. Consequently, we have 



H 



2m c H 

P 2 



P 2 C F a s C F p a 2 s (-§ + 2r/x - ±rV) 



2tt 



2m' + 



a 



+ Ar + Vl 5 + Vr + V hf + V rehK 



2m' c + + y L5 + V T + V h/ + Ke*,x, 



+ V LS + V T + V hf (3.13) 

(3.14) 
(3.15) 
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where 



H (o) 
A 



+ Ar . 



m c r 



4tt 



o 



Cfoi s 



_ 3 C F /3 a 2 s 
~ 2 2tt 

2n ' 



(3.16) 
(3.17) 
(3.18) 
(3.19) 



In the spherical coordinate, it is known that 

(i d 2 



H) 2 v 2 
i d 2 



-r H 

r dr 2 r 2 sin 9 89 



8 . a d 
smd de + 



o 2 



r 2 sin ( 



1 



XT 



18 8 Id 2 
sm0— + 



P? + 



r dr 2 r 2 \ sin 9 89 

L 2 

r 

We substitute Eq. ([330]) into Eq. ([3T6]) . and obtain 

1 

m, 



56> sin i 



(3.20) 



-— r + — ( J-^- s in0— + 1 02 

r dr 2 r 2 sin # 9$ 89 sin <9</) 2 



o 



+ Ar 



1 



L 2 



a 



— IP? + -xr - - + Ar . 



(3.21) 



A quantum state, with specified angular momentum quantum numbers I and m, is satisfied by 
L 2 \l,m) =l(l + l)\l,m) , (n\l,m) = Y l m {9,(t>), (3.22) 



where Y, 



are spherical harmonics. Therefore, in the calculation one can simply replace 



L 2 by its eigenvalue 1(1 + 1) in Eq. (|3.2ip . so that we have 



H (o) = ± 



P 2 , *(* + i) 



1 d 2 Z(Z + 1) 
xr + 7+ 



a 

- + Ar 

r 

a 

h Ar . 

r 



(3-23) 

m c r 8r A m c r A r 

Finally, we express H = H — 2m' c and concentrate on solving its eigenenergies. To perform the 
Pade approximation study, we decompose the Hamiltonian H = H — 2m' c into two parts, Hq 
and Hl- 



H = H-2m'=H c + H 



where 



H c = +V LS + V T + V hf 



(3.24) 
(3.25) 
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1 d 2 1(1 + 1) 

rr + 



(1-6) 



m c r dr 2 m c r 2 

1 d 2 1(1 + 1) 



^ + V LS + V T + V h f, (3.26) 



H L = + V rel>K (3.27) 



7 r + 



+ Ar + V reltK , (3.28) 



m c r <9r 2 m c r 2 

with < < 1. Here He is the Hamiltonian contains the Coulomb potential, while Hl involves 
the linear potential. 

B. Masses of charmonia obtained from the modified Pade approximation 

As the Hamiltonian exhibited in Eqs. f)3. 13|) and (|3,14p . at large distance, r > 1 fm, the strong 
interaction potential is expected to rise linearly, so that we cannot treat the linear potential as 
perturbative term to obtain the solutions. In the present study, we have one numerical parameter 
b, and three physical parameters: the charm quark mass m c , strong coupling constant a s , and 
the scale [i. We have checked that for n < 6 the intrinsic numerical errors are less than 7% if 
using 0.1 < b < 0.6. The intrinsic error measures the difference between the exact value and 
its Pade approximant. Putting the constraint on b, so that the intrinsic numerical errors are 
less than 2% for all states with n < 6, we use the method of the modified Pade approximation 
to approximate the masses of charmonia and then fit them with the current mass data for 
the S*-wave charmonium states, 77 c (l 1 5o), J/ip(l 3 Si), t] c (2 1 Sq), and ip(2 3 S\) to determine the 
remaining three physical parameters. After that we can determined the 5"-wave mass spectrum. 

To obtain the eigenenergies of H, we first define 

H(P) = H C + PH L , (3.29) 

where f3 is a real positive number. We can perturbatively solve the eigenenergies of H(f3) in the 
two limits (3 <C 1 and ^ > 1. As f3 = 1, the eigenenergies of H(f3 = 1) correspond to the real 
cc system. Once we have the results corresponding to (3 <C 1 and /3 > 1, we can interpolate 
the two limits to obtain the eigenenergies of the real cc bound states by using the method of 
the modified Pade approximation. Note that V r ei,K is contained in Hl and its contribution is 
well under control in the perturbative calculation in the limit /3 <C 1. In the calculation of the 
large /3 limit, V re ^K is perturbatively small compared to H L °\ If V re i t x is involved in He, its 
perturbatively correction to Hq will be 0(a 2 /(4b 3 )) and out of control for b ~ 0.275 which is 
obtained in the later study. For the S-wave charmonium system, the mass spectrum is described 
by the eigenenergies 2m' c + £^qq][2/1](1), where £'| 1 qq][2, 1](1) are the modified Pade solutions 
of H(f3 = 1) which will be explained below. 

For a small /3, using the Rayleigh-Schrodinger perturbation theory as given in the quantum 
mechanics textbooks 35], we obtain, for the S-wave states, 

E n00sW) = E nOol + E nollfi + ■■■ , (3-30) 



10 



where -E^ooi = -^nOOs is given by Eq. (|A24|) for which the detailed calculation can be found in 
Appendix O and 

-^nOOs = (V'nOOsl^lV'nOOs) 

tp^ 00s (r)H L ip n00s (r)r 2 sin 6drd6d<j) 



1 — b - 1 d 2 \ 

-R ns (r)-^rR ns (r)r 2 + XR ns (r)rR ns (r)r 2 ) dr . 



(3.31) 



Here R ns {f) is defined by Eq. (|A26p . and IV'Sms) are the eigenkets of the Hamiltonian He, 
given by 



l^nlms) = \tPnlms) ® \s, 8 Z ; S lt S 2 ) , 



(3.32) 



where \ipnims) and |s, s z ; s\, S2) respectively correspond to the spatial and spin parts of the wave 
functions and we have 



(r\lpn00s) = 1pn00s(r) , 

(see also Eq. (|A26|) for the detailed expression). 

For a large /3, we can rewrite the Hamiltonian in the following form 



(3.33) 



3(13) = p 



+ V rel , K + -{HP + V LS + V T + V hf ) 



(3.34) 



For the 5-wave states, the energy spectrum E> 00s corresponding to the large j3 limit can be 
expanded in power series with respect to 1//3: 



E„ 00.- < >') - P ( £^00s + ^nObs + ^nOCli^ + 



(3.35) 



Here i^jos i s the zeroth corrections, and -E^joi an< ^ ^nooi are ^ ne ^ rs * corrections. They can be 
evaluated by the perturbation theory and are 

"(l-^A 2 ' 1 



E. 



E 



>(0) 
nOOs 

>(1) 
nOOs 



(3.36) 



(0) W V re i, K ^ (0) (r)r 2 sin 0drd0# 



B„ 



At((£fcA)i 



3r + r n 



-1 a 4 



4 m 3 r ^ r 4 



Ai((^A)i 



.4vrr 2 (ir, (3.37) 



>(i) 

nOOs 



(^ (0) l^c|^ (0) ) 

^°\r)(HP + V h f)il>W (r)r 2 sin 6drd0dxt> 
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B, 



Ai((^A)i 



-6 d 2 



m c r dr 2 



rB r , 



B n 



Ai U^_X)zr + r n 



r 3m 2 



s(s + l) 



3\ 6(r) 



2 J 4vrr 2 



>4nr 2 dr, (3.38) 



where the contributions due to Vls and Vr vanish (see Appendix[Aj, are the eigenenergies 

of Hl with I = 0, and the corresponding eigenfunctions are 



^ (0) (r) 



3r + r n 



with being the so-called Airy function, the normalization: 



B r . 



[Mr]] 

(t^a) 



-1/2 



—Andr 



(3.39) 



(3.40) 



and r n the roots of the Airy function. 

Now we compute the energy spectrum for S-wave cc bound states using the modified Pade 
approximation. The eigenenergies of the Hamiltonian H(f3) are approximately by the modified 
Pade approximants: 

j&gp/n<ffl = m+ r?^ < 3 - 4i » 

The coefficients Po,Pi,P2, and gi can be determined in the following way. Comparing with 
Eq. (|3.30p . for a small (3, we have 

=Po~ (Pi -P0Qi)/3 + 0(f3 2 ) (3.42) 

= E<Sl + E<g>/3 + 0(/3 2 ) , (3.43) 

On the other hand, comparing with Eq. (I3.34H . for a large f3, we get 



P2 



P2.Pl 
1 ' 

Qi qi 

>(i) 



We therefore obtain the relations: 



P0 - AiOOs > 
Pi - PO^l = £ n00fl i 



1 





P2 



91 
Pi 



>(0) 
nOOs ' 



P2 PI _ „>(1) 
" ,2 ~ l ~ ~ ~~ -^raOOs 



(3.44) 
(3.45) 

(3.46) 
(3.47) 
(3.48) 

(3.49) 
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and arrive at the eigenenergies of real S-wave cc bound states: 



Sft2[2/l]03 = l) = ^±ZL±^, (3.50) 

J- + qi 



with 

,<(o) 



Po = Kooi , (3-51) 

<(0) >(0) _ p <(l) p>(l) 
- p<(0) _ „>(!) ' ^.52) 

P >(0)/ p <(l) , p>(0)\ 
P2 ~ <(0) _ „>(!) ' {0.06) 

_F<W 1 F>(°) 

91 " p<(0)_ F >(l) • ^.54) 

IV. NUMERICAL ANALYSIS AND DISCUSSIONS 

We have only three parameters, m c , a s , and fi to be related to the physical masses, and one 
numerical parameter b. The b is related to the intrinsic error for the Pade approximant compared 
to its true value. In the fit, we put the constraint on b, so that the intrinsic error of the modified 
Pade approach is small enough < 2% for the states with the radial quantum number n < 6. 
In general, the error is less than 7% for 0.1 < b < 0.6. (We will further discuss the intrinsic 
error later.) We adopt the masses of four well-measured 5-wave charmonium states, ?7 c (l 1 5o), 



J/tp(l 3 Si), n c (2 1 So), and ^(2 3 Si) 36(], as inputs to determine m c , a s , and \i. In Particle Data 
Group (PDG) [36f] , ^(2 3 5i) is denoted by ip(2S) or -0(3683), which will be discussed later. 



Under the intrinsic error < 2%, we perform the best x 2 fit which is defined by minimizing 

X 2 = xlo + X11 + X20 + X21 , (4.1) 

with 



X 



(expt t}i \ 

m n00s ~ rn n00s 
r expt I 
0m n00s I 



(4-2) 



where m e ^ s ±bm^ s are the experimental charmonium masses and m^gos = ^ m ' c Jl ~ ^nool [2/l](l) 
the theoretical predictions that we have calculated in this paper. We find that the minimum of 
X 2 is Xmin = 0-26 as well as b ~ 0.275 corresponding to the almost smallest intrinsic error in 
the fit. Our results, in good agreement with the data, can successfully account for the hyperfine 
splitting for the IS state as well as for the 25 state. The best fit values for parameters are 

mjji) = 1.517 GeV, (4.3) 
a 8 (n) = 0.273, (4.4) 
H = 1.843 GeV . (4.5) 
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One that note that m c , a s , and fi are not really physical parameters since an isolated charm quark 
cannot be observed. The renormalization scale fj,, which is adopted to separate the potential 
into several parts, is chosen to be the quantity that charmonium becomes stable, so that after 
some combination we have m c — > m' c and a s — > a s as shown in Eqs. (3.12) and (3.13). 

It is interesting to note that the Coulombic and linear potentials defined in Eq. (|3.14p are 
then obtained to be 

0.288 



a 

h Xr 

r 



+ (0.241 GeV^) r, 



(4.6) 



while in the Cornell potential model the potential is parametrized as V = —a/r + er, with 
a ~ 0.52 and e ~ 0.18 GeV 2 [10]. For the Cornell potential, a is usually identified by CfO- s . 
Nevertheless, our a is 



C F 



1 , ai+7£/3 /2 
1 H a s 



7T 



3 C F f3 a 2 s 
2 2vr 



Using the obtained parameters, we can further get the masses of higher S"-wave states. The 
results are given in Table HH For c omp arison , we also list the current data assignments j^] and 
some other theoretical results 13l. Il4. 137H39I]. It was known that ■0(3683) and ^(3770) could 



be the mixtures of 2 3 Si and l 3 Z?i states due to the fact that, instead of L, the total angular 
momentum J is a conserved quantum number; L can be broken by some relativistic effects, for 
which especially the operator of the tensor force does not commute with I? for states with S = 1 
(see the discussions in Ref. [2l|). However, we see that the fitted Xmin ^ s consistent with zero, 
which may hint that the S-D mixing effect is negligible. One should note that the relatively large 
e + e~ width of the ip{2>77{)) is difficult to understand if it is a pure 1^D\ cc state. This problem 
can be solved if "0(3770) has an admixture of 2 3 5*i cc state 15]. Our results show that if the 
smallness of S-D mixing effects can be applied for higher radial excited states, the singlet-triplet 
splitting mass difference is about 50 MeV for n > 3. Conventionally, the observables ^(4040) and 
^(4160) were assigned as 3 3 Si and 2 3 Z?i states, respectively. However, our calculation suggests 
that the ^(4160) may be dominated by the 3 3 Si state. It has been noted that the ^(4160) has 
a much larger e + e~ width so that it may have a significant S- wave cc component 40]. In the 
flux-tube model, the light hybrid charmonium states lie ~ 4.1 GeV and it was suggested that 
■0(4040) and ^(4160) may be the strong mixtures of the hybrid charmonium and i>(3S) 411 ] . 
which can explain why T e+e ~ (^(4040)) ~ T e+e ~ (^(4160)). If so, the X(4160) might be further 
assigned as the r] c {3 1 S 1 ). The X(4160) with a mass of (4156 ± 29) MeV/c 2 and a total width of 
r = (139^g5 3 ) MeV/c 2 was seen by Belle in the DD* recoiling from the J ft/) in the annihilation 
process e + e~ — > J/tpD*D* 42|. We obtain the mass for the 4 3 Si to be 4475.1 ±13.5 MeV, which 
is about 50 MeV larger than ^(4415) which is conventionally assigned as the 4 3 Si state. These 
discrepancies can be further clarified by including higher order corrections in the calculation. 
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TABLE II: The predictions for the S-wave charmonium spectroscopy (in units of MeV), com- 
pared with the results from the current data assignments (PDG) and from theoretical calcula- 
tions by lattice QCD (Lattice), perturbative QCD-based (PQCD), and phenomenological quark 
model (QM). The theoretical errors are estimated in Table Unl 



State (n zs+1 Lj) 


PDG 36J 


ml ■ 1 

1ms work 


T j_ i ■ Too onl 

Lattice [38, 39J 


PQCDI37I 


QM|13J 


QM[14] 


_ /lie \ 


2980.5 ± 1.2 


onon c inn 

2980.5 ± 0.0 


om A 1 A 

3014 ± 4 


on cr a 

305b 


2979 


onnn 

3000 


J/^(1 3 S 1 ) 


3096.916 ± 0.011 


3096.9 ±0.0 


3084 ± 4 


3097 


3096 


3100 




3637 ± 4 


3634.9 ± 3.0 


3707 ± 20 




3583 


3670 


^(2 3 ^) 


3686.093 ± 0.034 


3686.1 ±3.0 


3780 ± 43 




3686 


3730 


Vci^So) 




4068.4 ± 8.0 






3991 


4130 


</>(3 3 Si) 


4039 ± l a 


4118.6 ±8.0 






4088 


4180 


r/ c (4 1 5 ) 




4424.5 ± 13.5 










V>(4 3 Si) 


4421 ± 4 & 


4475.1 ± 13.5 








4560 






4730.0 ± 18.8 










</>(5 3 Si) 




4781.8 ± 18.8 














5000.4 ± 23.9 










<K6 3 Si) 




5052.0 ± 23.9 











a It is called t/>(4040). 
6 It is called V(4415). 



For estimating the numerical uncertainties for our predictions , we take into account the 
eigen-solutions for the spatial Hamiltonian H^> with I = 0, i.e., the full Hamiltonian H without 
the spin-spin interaction term and with 1 = 0: 

(- — f- 2 r ~- + Xr) *$> = £<M°> • (4.7) 
V m c r or z r I 

As in Eqs. (|3.26p and (|3.28|) . we introduce the parameter b and then split into two parts, 
Hq^ and Hl. Substituting the values for m c , A, a, given in Eqs. f)4.3[) . (|4.5p . (|4.6p and adopting 
b = 0.275, we numerically solve the above equation. Table IIIII compares the modified Pade 



Without existence of the spin-spin (hyperfine) interaction, one could solve numerically the Schroedinger 
equation for the "spatial part" of the Hamiltonian H, i.e., determine its spatial wave functions. How- 
ever, it is highly nontrivial if there exists the spin-spin interaction and one would like to fit numerical 
solutions to the data. The spin-spin interaction is relevant to explain the mass splitting between 
?7 c (l 1 S , o) and J/tp(l 3 Si) and between ■q c (2 1 So) and J/ip(2 3 Si) > which cannot be computed well in the 
literature so far. 
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results with numerically exact eigenenergies. We see that the modified Pade approach yields 
approximations which are not larger than 1.5% of the exact solutions for states with n < 6. 

TABLE III: Error estimate in the modified Pade approach. The eigenenergies are in units of 
MeV. 



State (nS) 


e£ } (Exact) 


E ( n } (Pade) 


Pade-Exact (error) 


(Pade-Exact) /Exact 


IS 


603.26 


603.29 


0.03 


0.01% 


2S 


1253.87 


1256.86 


2.99 


0.2% 


3S 


1761.25 


1769.20 


7.95 


0.5% 


4S 


2202.15 


2215.61 


13.46 


0.6% 


5S 


2601.77 


2620.58 


18.81 


0.7% 


6S 


2972.23 


2996.09 


23.86 


0.8% 



V. SUMMARY 

The 5-wave charmonium spectroscopy has been calculated by considering the Hamiltonian 
with the non-relativistic QCD potential. For the next-to-leading order QCD loop corrections to 
the potential, we expand the logarithmic factor In \ir about r = l//x, where [i corresponds to 
the typical charmonium scale of order m c v, so that the QCD potential can be modeled as the 
Coulomb plus linear form, which is consistent with the Cornell potential. In our approach, we 
have performed the best x 2 fit by comparing the current mass data of the S'-wave charmonium 
states, r/ c (l 1 5o), J/ip(l 3 Si), / q c (2 1 So), and ip(2 3 Si), with their modified Pade approximants. Our 
results, in good agreement with the data, can successfully account for the hyperfine splitting for 
the 15 state as well as for the 2S state. The fitted parameters are m c (fi) = 1.517 GeV,a s (fi) = 
0.273, and jj, = 1.818 GeV, consistent well with the ranges that one usually used. Using then 
these three parameters we have further predicted the S'-wave mass spectrum with n < 6. 
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Appendix A: Eigenenergies and eigenfunctions of the Hamiltonian He 



The eigenenergies and corresponding eigenfunctions of the Hamiltonian operator He can be 
calculated by using the Rayleigh-Schrodinger perturbation theory. We decompose He into two 
Hermitian parts, H^P and the rest, 

H C = H§ ) + (V LC + V T + V hf ), 



where 



V LS 
V T 
Vhf 



3C F a s 
2 m 2 r 3 

C F a s 
4 m 2 r 3 



L-S. 



'12 



C F a s f ^f-S^-f) - - N 

o 9 Ol • £>2 



4 m 2 r 3 



8irC F a s 4nC F a s 
-Si ■ S 2 5°(r) = 



3m 2 : 



3m? 



8(8 + 1) 



5 3 (f) 



(Al) 

(A2) 
(A3) 
(A4) 



To solve the eigenvalue problem of He, we instead consider the following Hamiltonian function 

H c (k) = H { £ ) +K (y LC + V T + V hf ), (A5) 
so that we have He = He(l), where k is a continuous real parameter. 



1. The corrections due to the spin-obit and tensor interactions in S*-wave states 
of Hamiltonian He 



It is obviously that (Vls) = 0, where (• • •) means the expectation value of the S- wave states, 
while (Vt) is given by 

C F a s 



(V T ) 



4 m 2 r 3 



(S 



12 



C F a s ^ (r • Si)(S 2 -f) _ ^ _ ^ 



4 m 2 r 3 



(A6) 



where 



(3 (r 1 5 1 )0Vr2 _ ^ ^ 



= (J d 3 r\^ n00s (r)\< 
= (J d\\^ nms (r)f 



] 3(r-Si)(S 2 -r)-S l -S 2 r 2 



3(Si x x + Siy y + S lz z)(S 2x x + S 2y y + S 2z z) Si ■ S 2 r 2 



KA7) 
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with V'nOOs(r) being the spatial part of the wave function for the Hamiltonian He, given in 
Eq. (|A26|) . which is spin-dependent but independent of (f> since m = 0. In terms of the variables 
of spherical polar coordinates, r, 6 and (j>: 

x = r sin 9 cos <f> , (A8) 
y = r sin 6 sin <f> , (A9) 
z = rcos9, (A10) 

we have 

y xydfi = J r 2 sin 3 cos <f> sin (f>d9d(j) 

= = J yzdQ = J xzdn. (All) 
x z dQ. = I r 2 sin 3 6* cos 2 



A jr 2 = I y 2 dfl = I z 2 dU. (A 12) 



Using the results of Eqs. (|A11|) and ()A12p to the angular integral in Eq. (|A7p . we obtain 
J dQ(r- S±)(S 2 ■ r) = J dU(Si x x + S ly y + Si z z){S 2x x + S 2y y + S 2z z) 
= J dQ-r 2 (Si x S 2x + Si y S 2y + Si z S 2z ) 



dn-r 2 S! ■ S 2 , (A13) 
3 



and therefore 

(Vt) 



C F a s / 3 2 3(S lx x + Si y y + S lz z)(S 2x x + S 2y y + S 2z z) - Si ■ S 2 r 2 

CpOs f , 3 . i2 £h ■ S 2 r 2 - Si ■ S 2 r 2 

J d r \^Ws\ ~ 5 ) = 0. (A14) 



2. The corrections in S-wave states of Hamiltonian He 

Since the spin-orbit and tensor interactions can be neglected for the S-wave states, following 
the standard approach, the eigenenergies E n i ms (n) of Hc(k) can be determined in terms of the 
perturbation expansion: 

E nlms {K) = E^ s + eQ^k + E$ ms K 2 + 4/ + • • • , (A15) 
(0) _ _nW±_ , 
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E 



(i) 

nlms 



C (o) MC F a s 

\ ^nlms I 



3m 2 . 



sfs + 1) 



-J 

Jo 



VramW 3m 2 

3m? 



a(a + l)- 



«(s + l)-2 



<5 3 (»=)^Sm ( f ) r2 sin Odrd0d(f> 
(2l + l)^^R^{r)r 2 dr, for m = J , 



(A17) 



(2) 
nlms 



Vnklms'knlr 



h(0) 
J nlms 



klms 



E 



(3) 
nlms 



E 

fc=i,... 

E S , p (0) _ „(0) ,,,10) _,,(()) ^— ,,.,(()) m, 

fc=l,... g=l,... K-^nlms ^ klms iK 1 ^ nlms 11 'glms> *=}•- K^nlms ^klms) 



^nklms Vkglms Vgnlms 



(A18) 



~77M Tn(o) ^ 2 A Aiy J 



where E^ ms and IVv^ms) are respectively the eigenenergies and eigenkets of -f^, and 

Vnklms = V kn i ms 

= ^nlms\ n o Si ■ S 2 5 (^fclLi) 



3m^ 



The states \^Zms) read 



«(« + l)"2 



(2/ + l)^^R { ^(r)r 2 dr, for m = ].(A20) 



l^nlms) = I^Sm)®l s ' s *; a l' a 2), 



(A21) 



where 



iT(M) 



nb 



1 (n — Z — 1)! / m c ar\ am c 2 i+i 

1 g 2nb ±j 



n-l-l 



m c ar 
nb 



(A22) 



V 2n[(n + /)!] V n& 

In the calculation, we introduce the transformation, called the Pade approximation, to accelerate 
the convergence of the perturbative series of E n i ms (n) which are approximately presented as 
rational functions, £^^^ s [2/1](k), 



E n lm S (K) ~£^[2/1](k 



PO 7 +Pi'k+P2 , « 2 



(A23) 



n(3) 



1 + 9p/e 

where pg 1>2 and gf can be determined by the values of E$ ms , E$ ms , E$ ms , and E^' ms . 

Therefore, for the 5-wave states, the eigenenergies of the Hamiltonian He are approximately 
to be 

nP J_ JO _j_ JO 

(A24) 



Enoos — E n00s (l) ~ S^oo [2/l](l) — ^° : ^ c P2 , 



and the corresponding wave functions are 

IV'nOOs) = IV'raOOs) ® |S,S*;S1,S2) 



(A25) 
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where up to the second order 

/T"r„, m , x 



with 



47T 



'nOOs v) 

4 )(r) + i?W(r) + 4 2 s )(r); 



= E 



ma 3 /(n + 1)! fmar\ am c . r t fmar s 
~^h'\ 2n(n!) 6 ""^ n_1 \~n&~, 

-^ C(0) f 



^ p(0) _ p (o) %0 ° 1 j 

fc=i,... -C/ n0 0s -^fcOOs 



^SSw = E E 

fc=l,... 9=1,... 



Kfcoo, „(o) , v fT _ /T (i) , , 

^ p(0) _ p (0) ^ 1 j V4vr ~ V47T nsl j ' 

fc=i,... £/ n QOs -^fcOOs 

V"^ VnaaOxVakOOis ,0(0), 



k= 



(F (P) _ p(0) 



- E 

*:=i 

T 



ngOOsVgfcOOs , C(0) , -* 

j^fcOOs^nnOOg ,C7(Q),-* 
- / p (0) _ p(0) s 2 ^00 W 

1 



^flS(r). 



(A26) 
(A27) 



(A28) 



(A29) 
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